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Восстановление формы объемных 
медицинских объектов с помощью 

„Ж* 
дистанционных карт поперечных сечений 


Представлены новые методы для реконструкции трехмерной поверхности объекта из нескольких 
замкнутых, в общем случае неплоских кривых, включая контуры, которые были очерчены вручную. 
Для восстановления объектов разной формы используются двухмерные дистанционные карты. 


Введение 


Задача реконструкции трехмерной поверхности объекта хорошо известна в раз- 
личных областях, особенно в медицине. Она используется при работе с изображе- 
ниями магниторезонансной томографии (МРТ), компьютерной томографии (КТ), 
трехмерного ультразвукового исследования (30 УЗИ). 

Задачи реконструкции можно разделить по уровню взаимодействия на авто- 
матические, полуавтоматические и ручные. В общем случае задача может быть сфор- 
мулирована как восстановление объекта из произвольного облака точек. 

Но обычно работа осуществляется с двухмерными контурами, представля- 
ющими области интереса на сечениях медицинского изображения. Получение этих 
сечений может быть сложной проблемой, и оно может производиться вручную или 
автоматически. Например, оконтуривание органа на ЗО УЗИ - очень сложная проб- 
лема, и эта задача обычно выполняется медицинским экспертом. 

Существуют также некоторые особенности, связанные с размещениями контура 
в пространстве. Например, на МРТ и КТ контуры расположены в строго параллель- 
ных областях, а для некоторых типов УЗИ слои в зависимости от протокола скани- 
рования находятся в произвольных плоскостях. В некоторых случаях настоятельно 
рекомендуется восстановить объект лишь по нескольким слоям, особенно в случае 
ручной обводки. Это избавит медицинских работников от рутинной работы выделе- 
ния каждого слоя. 

В последние годы предложено много методов реконструкции. Прежде всего, стоит 
упомянуть растровую реконструкцию объектов. Эти методы обсуждены в работах 
Г. Трис [1], [2]. Статьи описывают различные алгоритмы: для параллельных слоев [2] 
рассмотрены случаи объектов сложной формы. 

Основная идея этой группы методов — особый подход для межслойной вок- 
сельной интерполяции, зависящей от взаимного расположения слоев, а также формы 
поперечного сечения. 

Нужно отметить, что такие воксельные методы обычно медленны и неудобны 
для использования и визуализации. По этой причине разработано много методов для 
аппроксимации поверхностей объектов сплайнами и кривыми различных типов. 


у Работы выполнялись в рамках проекта [ТС В-1486. 


370 «Искусственный интеллект» 372010 


Восстановление формы объемных медицинских объектов... дн 


Например, можно восстановить векторную поверхность объекта неравномерными 
рациональными би-сплайнами [3], которые хорошо аппроксимируют обведенные кон- 
туры объекта. Но качество этого подхода сильно зависит от числа контуров. И не- 
всегда возможно восстановить контуры объекта для всех слоев, особенно в случае 
ручной обводки. 

В случае объектной реконструкции по нескольким контурам можно исполь- 
зовать один из методов на основе модельных представлений, которые работают 
лучше на определенных типах объектов. Один из таких подходов использует сплайны с 
одним параметром для реконструкции щитовидной железы на трехмерных УЗИ [4]. 

В другой работе [5] рассмотрена более общая проблема — реконструкции формы 
от произвольных поперечных сечений. Предложенный подход состоит из двух ос- 
новных шагов: вычисление расположения секущих плоскостей и реконструкция, и 
аппроксимация объекта из его пересечений с границей ячейки. Возможная проблема 
ветвления решена использованием триангуляции Делоне. Этот метод довольно слож- 
ный, однако он может быть применен в приложениях реконструкции. 

Например, чтобы восстановить простой медицинский трехмерный объект из непарал- 
лельных поперечных срезов Богуш и др. [6] используют двухмерный морфологический 
морфинг, аффинные преобразования и интерполяцию сплайнами третьего порядка. 

Также существует большой класс методов для восстановления формы объекта 
из произвольного облака точек. Этот класс методов используется в различных задачах, 
таких, как дистанционное зондирование, стереоприложения и т.д. Но эти подходы 
также могут быть использованы в медицинских приложениях, особенно на входных 
данных высокого разрешения. 

Большой подкласс таких методов — алгоритмы симплекс-сеток. В [7] представлены 
методы симплекс-сеток для восстановления трехмерных ультразвуковых изображений. 
Алгоритмы используют непараллельные сечения для симплексной адаптации сеток 
из вручную обведенных границ объекта на нескольких представителях. Подобный 
метод представлен в [8]. 

Также интересны следующие способы восстановления облака пикселей — алгоритм 
катящегося шара [9] и байесовская реконструкция поверхности [10]. Первый из них 
является своего рода морфологическим восстановлением на векторном пространстве. 
Второй — использование статистического метода для реконструкции и последующей 
фильтрации шумов на поверхности. 

Постановка задачи. Наш метод используется для тех же задач и основывается 
на формировании промежуточных слоев на основе дистанционных преобразований. 
Он достаточно прост в реализации и понимании, и работает для сложных медицин- 
ских изображений. 

Главное преимущество метода — возможность генерировать произвольное коли- 
чество промежуточных слоев. Кроме того, предложен оригинальный подход к решению 
проблемы ветвления. В общем случае, алгоритм может работать на параллельных и 
непараллельных слоях. 


Определение контуров объектов для промежуточных слоев 


Построение промежуточных слоев может быть реализовано при помощи ана- 
лиза свойств дистанционных карт [11-13]. Если число объектов начального слоя такое 
же, как и конечного, то построение промежуточного слоя довольно простое. Но в 
случае ветвления эти методы имеют ряд проблем. 

Тем не менее, свойства дистанционных карт могут быть использованы для пост- 
роения промежуточного слоя. Но в этом случае будут построены два промежуточ- 
ных слоя. Рассмотрим случай, когда первый слой содержит один «корневой» контур, 
а второй слой — несколько контуров (рис. 1). 
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Рисунок 1 — Последовательные бинарные слои: а) с одним целым объектом, 
6) с разделенными объектами 


Необходимо использовать одинаковое количество контуров для начального и 
конечного слоя. Для решения этой проблемы используются алгоритмы построения 
выпуклой оболочки, однако все они работают только для связных объектов. Таким 
образом, необходимо соединить все объекты на слое. Объединение объектов реали- 
зовано через вычисление координат центра масс и проведения линий через них. 
Затем для всех объектов строится выпуклая оболочка (рис. 2). 


Рисунок 2 — Объединение объектов: а) соединение центров масс, 6) построение 
выпуклой оболочки 


Область определения контура объекта для среднего слоя заключена между грани- 
цей начального слоя и выпуклым контуром конечного слоя (рис. За). Эта область стро- 
ится путем вычитания слоев. Этот контур находится из дистанционной карты (рис. 36), по 
которой можно построить равноудаленный от границ замкнутый контур (рис. Зв). 


а) 6) в) 
Рисунок 3 — Определение контура первого промежуточного слоя: а) область для 
определения контура промежуточного слоя, 6) его дистанционная карта и в) линия раздела 


Следующая задача — разделить объект на среднем слое на несколько. При этом 
число объектов соответствует числу на конечном слое. Линии разделения должны 
быть проанализированы с позиций объектов на конечном слое. Лучшее решение, 
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когда они используют линии водораздела между объектами конечного слоя (рис. 4а). 
Вычитание этого водораздела позволит построить второй промежуточный уровень 
(рис. 46). Этот слой включает то же количество объектов, что и на конечном уровне. 

Следующие шаги будут по определению слоев от начального до первого сред- 
него слоя и от второго среднего слоя до конечных слоев. Обработка слоев идентична, но 
второй случай более сложный. Поэтому при следующем описании мы будем использо- 
вать только построение на этапе между вторым средним слоем и готовым слоем. Для 
этого этапа каждый слой содержит одинаковое количество контуров. Однако объекты 
имеют различный размер, положение и форму. Задача состоит в нахождении оптималь- 
ных контуров объектов для промежуточного слоя. Эти контуры соответствуют оптималь- 
ному расстоянию между первым и вторым слоями. Это легко определить, используя 
конъюнкцию и дизъюнкцию изображений слоев (рис. 5). 


_ 


а) 6) 


Рисунок 4 — а) линии водораздела для последнего слоя и 6) объекты на втором 
промежуточном слое 


6) 


Рисунок 5 — а а) конъюнкции и 6) дизъюнкции второго 
промежуточного и последнего слоев 


Разница полученных результатов позволяет найти область определения контура 
(рис. ба). Конъюнкция этого результата с контуром второго промежуточного слоя (рис. 66) 
позволяет сохранить связность для построенных контуров. 


а) 6) 


Рисунок 6 — а) область определения контура второго промежуточного и последнего 
слоев, 6) контуры второго промежуточного слоя 
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Следующий алгоритм (рис. 7) позволяет определить объекты с контуром области 
и несвязанные объекты. 
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Рисунок 7 — Алгоритм классификации объектов для построения контура 


Проблема заключается в определении областей при построении дистанционных 
карт в двухмерном пространстве. На изображении объединенных слоев возможно 
как связанные, так и разделенные объекты. Области определения контура строятся с 
помощью алгоритма, блок-схема которого представлена на рис. 8. 

Для всех разделенных объектов анализируется расстояние между центрами масс. 
На основе этой информации формируются пары объектов. Для каждой пары строится 
выпуклая оболочка (рис. 96), которая затем корректируется с помощью исходных 
данных предыдущего уровня (рис. 9в). Объединение исправленного выпуклого контура 
и объектов позволяют определить области для определения контуров. 
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Рисунок 8 — Алгоритм построения области определения контуров 


а). 6) 


в) 


Рисунок 9 — Обработка разделенных объектов: а) разделенные объекты, 6) выпуклая 
оболочка, в) коррекция границ, г) новый связанный объект 


Дизъюнкция результатов обработки связанных и несвязанных объектов позво- 
ляет определить промежуточные контуры с помощью водораздела. Аналогично опреде- 
ляются контуры объектов для остальных промежуточных слоев. После создания 
набора промежуточных слоев возможно провести трехмерную реконструкцию с более 
высоким разрешением (рис. 10). 
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Рисунок 10 — а) области определения контуров для промежуточного уровня, 
6) объекты промежуточного уровня 


Результат работы алгоритма для входных данных на рис. 1 представлен на рис. 11. 


Рисунок 11 — Результат трехмерной реконструкции с помощью 20 промежуточных слоев 


Заключение 


Разработанный алгоритм позволяет получить качественный результат благода- 
ря восстановлению промежуточных сечений. Использование связывания всех объек- 
тов на слое обеспечило стабильность работы алгоритма для объектов любой формы, 
а растровое представление данных исключило элементы скручивания поверхности. 
Использование дистанционной карты расстояний позволило быстро и эффективно 
определить реальный контур между слоями и сохранить значения всех исходных рас- 
стояний. Поэтому данный алгоритм наиболее эффективен для использования его в 
задачах, включающих измерения объектов. 

Полученные результаты будут способствовать развитию диагностического прог- 
раммного обеспечения лучевых методов для правильного расчета данных по динамике 
регрессии опухолевых образований и размеров резидуальной опухоли, оптимизации при- 
нятия решения по вопросам лечения. Определение правильного пути лечения позволит 
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значительно снизить стоимость лекарств и лечебных процедур. Работа направлена на 
улучшение информационного содержания данных первичного диагноза, а также для обес- 
печения комфорта специалистам-медикам в ходе анализа результатов обследования. 
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О.М. Нядзьведь, П.В. Лукашев, Д.А. Гончаров 

Виновлення форми об’емних медичних об’ектв за допомогою дистанцйних 

карт поперечних перетинв 

Представлен! нов! методи для реконструкций тривимрно! поверхн! об’екту з деюлькох замкнутих, в 
загальному випадку неплоских кривих, включаючи контури, як! були обкреслен! уручну. Для вдновлення 
об'ект1в р1зно! форми використовуються двом1рн! дистанщйн! карти. 


А.М. Меасуеа, Р.Т. Гиказйеясй, Р.А. Напсйагои 

Весопзегисйоп о7 Мед са! ОБесЕ Зваре Ву О апсе Марз оЁ 5 Сго$$ ЗесНоп$ 

ТБе пеуу {есшаиез 15 ргезеге4 ш ог4ег ю гесопзиласЕ ЗО оБ]есё зигРасе от зеуега| с1озе4, т зепега|, 
поп-р]апаг сигуез, шсте сопюиг$ фай \уеге оц тед тапиаПу. 20 д1апсе тар изе4 №0 гесопзёласЕ оБуес{ оРа 
Чегет $Варез. ВтапсЬ те рго ет а150 91$са5$е4. 
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